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Abstract 

We review recent progress in applying relativistic hydrodynamics to the modeling of heavy-ion 
collisions at RHIC and LHC, with emphasis on anisotropic flow and flow fluctuations. 



1. Introduction 

Relativistic hydrodynamics is the only first-principles approach |[T1 to the out-of-equilibrium 
dynamics of the strongly-coupled quark-gluon plasma formed in heavy-ion collisions. It is a 
macroscopic description in which the plasma is modeled as a continuous lump of fluid expanding 
into the vacuum. It plays an important role in understanding the soft sector of nucleus-nucleus 
collisions at RHIC and LHC. 

The goal of relativistic hydrodynamics is to describe the bulk of particle production. A 
hydrodynamic calculation uses as input a model for initial conditions — typically, for the initial 
density profile in a collision. The only input needed for hydrodynamics is the initial value of 
the energy-momentum tensor: all other microscopic details of the initial state are washed out by 
local thermalization, which is implicitly assumed by the hydrodynamic description. One then 
solves for the fluid expansion using equations of ideal or viscous relativistic hydrodynamics. 
Eventually, the fluid is transformed into independent hadron^ 

The fluid is continuous: single-particle distributions from the fluid can be computed (as a 
function of transverse momentum /?,, pseudorapidity 77, azimuthal angle (p) with unlimited ac- 
curacy for a given initial condition, unlike in an actual experiment. As we shall see later, this 
peculiar feature of the fluid description is an important one. Many groups implement hadroniza- 
tion using a Monte-Carlo generator, thus mimicking the experimental situation where at most a 
few thousand particles are observed per evenj^ This amounts to picking randomly one element 
out of a thermal ensemble. The fluid is the thermal ensemble. 

In 2010, it was shown that observed azimuthal correlations between particles separated by a 
gap in pseudorapidity (usually referred to as "long-range" correlations) are compatible with this 
hydrodynamic picture ||2][3l, thus providing a unique signature of collective behavior. In these 
proceedings, we review the status of this "flow hypothesis" in light of recent measurements. We 
then hst open issues, and recent works since the last Quark Matter conference. 



Depending on the implementation, further hadronic decays and/or rescatterings may occur after hadronization. 
^This is the most natural way of implementing a "hadronic afterburner", i.e., further rescatterings after hadronization. 
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Figure 1 : (Color online) Left: initial energy density in the transverse plane at z = in a random central {h = 0) Pb-Pb 
collision at 2.76 TeV simulated by Nexus. Right: azimuthal distribution of outgoing particles after evolving this initial 
distribution through ideal hydrodynamics, in several bins of transverse momentum p, (in units of GeV/c). 



2. A close look at a hydro event 

Quantum fluctuations in the wavefunctions of colliding nuclei result in an initial profile which 
is not smooth and fluctuates from event to event |4|. In order to understand the physics associated 
with these initial fluctuations, it is instructive to look at a particular event in some detail. We 
simulate a central {b - 0) Au-Au collision using the event generator NeXus [5J, which gives the 
initial density profile shown in Fig. [T] (left). These initial conditions are then evolved through 
ideal hydrodynamics [6]. As explained above, thermalization thus transforms a single event 
into a thermal ensemble. The term "hydro event" usually refers to this thermal ensemble. We 
transform the fluid into hadrons using a Monte-Carlo generator. This hadronization is repeated a 
few thousand times so that we can reliably calculate ensemble-averaged quantities for the single 
hydro event displayed in Fig. [T] (left). 

A central collision between spherical nuclei is azimuthally symmetric, except for quantum 
fluctuations. One clearly sees in Fig. [T] (left) that fluctuations break rotational symmetry in the 
transverse plane. Fig. [T] (right) displays the azimuthal distribution of charged particles in the 
pseudorapidity interval It/I < 1 for various bins in p,. Anisotropics are at the % level at low p, 
and become stronger and stronger as p, increases. 

The azimuthal distribution can be expressed as a Fourier series: 

^ «=-oo 

where V„ - {e'"'^} is the nth (complex) Fourier coefficient, and curly brackets indicate an average 
over the probability density in a single event. Writing V,, - v„e'"*", where v„ is the (real) 
anisotropic flow coefficient and the corresponding phase, and using y_„ = V* (where the 
superscript * denotes the complex conjugate), this can be rewritten as 

^ ^ = 1 + 2 ^ v„(p,, jj) cos n (0 - "Vnipt, rf}) . (2) 

^ n=\ 

2 




Figure 2: (Color online) Anows representing the magnitude and direction of the first Fourier coefficients of the curves 
in Fig.[T](right). The length of the arrow is Vfj and its orientation is ^fj. From left to right: n — \ (directed flow), /t — 2 
(elliptic flow), ;i = 3 (triangular flow). 



Note that, for this form to describe an arbitrary distribution, both v„ and *F„ may depend on trans- 
verse momentum p, and pseudorapidity 77. Most models of initial conditions predict fluctuations 
in the form of longitudinally-extended flux tubes [7|. Therefore, one expects and v„ to depend 
little on rapidity. Although this rapidity dependence is worth investigating ll8ll9l [T0ll . we focus 
here on the stronger pt dependence. 

Note that event-by-event anisotropic flow is thus well defined in hydrodynamics only, be- 
cause thermalization transforms a single event into a thermal ensemble. In an actual experiment, 
the relative statistical error is typically 50% for V2 (in a mid-central Pb-Pb collision analyzed by 
CMS or ATLAS), and 100% for V3. Therefore the event-by-event v„, as defined above from the 
single-particle distribution, cannot be measured experimentally. 

Fig.[2]displays the magnitude (v„) and direction (*P„) of the first three harmonics obtained by 
aFourier decomposition of the distributions in Fig.[T] Note that *P„ is defined modulo 2nln, and is 
therefore represented using n arrows pointing into the corresponding directions. As noted before, 
azimuthal anisotropies generally increase with p,, hence v„ increases with p,. The directions of 
elliptic flow (*I'2) and triangular flow Q¥t,) depend slighly on pj^ The p, dependence is much 
stronger for . This reflects the fact that the net transverse momentum of the central rapidity 
slice is expected to be close to 0, hence the integral of p,V\ (p,) = p,vi (/5,)e'*' over all particles 
should be 0. One therefore typically expect that *Fi rotates by n as transverse momentum goes 
from to (xj3 

3. Flow in data: a close look at 2-particle correlations 

As explained above, event-by-event anisotropic flow is well defined in hydrodynamics only. 
The number of particles per event in an actual experiment is too small to measure V„(pi, 77), or 



^The particular event that we randomly picked for this illustration turns out to have a rather small i'2 and large V3. 
''The hydrodynamic code (NeXSPheRIO) that we use has a small residual net transverse momentum due to the choice 
of initial conditions, as noted in a previous study of directed flow ilii . This eft'ect has not been coiTected here. 
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even the integrated v„. Anisotropic flow can only be measured through event-averaged azimuthal 
correlations between particles. The simplest azimuthal correlation is the pair correlation. One 
forms all possible pairs of particles and then performs a harmonic decomposition of the distribu- 
tion of the relative angle: 

V„A(/,?,;,f) = <cos(«(0„ - cPt))) = <e'«(*«-^">>. (3) 

This quantity can be measured as a function of the transverse momenta of both particles, labeled 
as a and b, thus yielding a correlation matrix. This correlation matrix has recently been measured 
by LHC experiments lfT2l [131 IT4 1 . It is symmetric by constructiorj^ 

We now evaluate this quantity in hydrodynamics and show that it is can be expressed sim- 
ply in terms of anisotropic flow. In a single hydro event, particles are emitted independently. 
Therefore 

^ } j^-™/} ^ yayh* ^ vyle"<r.-^''\ (4) 

where { ■ ■ ■ ) denotes an average over a single hydro event. The first equality in Eq. Q expresses 
mathematically that particles are independent: this implies that the two-particle correlation, when 
written in complex form, factorizes in a single hydro event. 

The experimental quantity, Eq. ([3]l, is then obtained by averaging over a larger number of 
hydro events: 

v„A(p^p?) = (Kv^) = (v>;y"<*"-*"') (5) 

Due to parity symmetry, only the real part remains after this average, hence the cosine in Eq. ([3]l. 

From this relation alone, one can make the following general statements about the event- 
averaged correlation matrix when flow is dominant: the diagonal elements must be positive, and 
the off-diagonal elements must satisfy a Cauchy-Schwarz inequality L15i . 

ynKip^p") > 0, 

\VnK{p%p'f < V,Ap",p';)V„a(p';,p'^). (6) 

It is often stated that flow implies factorization. As shown above, factorization holds for a single 
hydro event. It implies that the second inequality in (|6]l is saturated, i.e., equality is achieved. 
After averaging over events, one does not in general expect factorization to hold, therefore the 
2nd inequality is generally a strict one for p" + . Any violation of (j6| is an unambiguous 
indication of the presence of non-flow correlations. 

Inequalities (j6]) can be directly tested on experimental data. We use ALICE data |[T2l for 
0-10% central Pb-Pb collisions. For n = 2 (Fig. [3] middle), all diagonal elements are posi- 
tive, and the Cauchy-Schwarz is verified except for two matrix elements, where the violation is 
compatible with a statistical fluctuation (the deviation is barely above 1 cr). The ALICE col- 
laboration concluded from their analysis that the correlation matrix factorizes approximately for 
trigger particles below 4 GeV/c. However, deviations from factorization are clearly seen at much 
lower transverse momentum. It is interesting to note that they are everywhere compatible with 
the general inequalities predicted by hydrodynamics. If hydrodynamics is a valid approach, such 
deviations are a direct evidence of event-by-event fluctuations. 

For n = 3 (Fig. [3] right), diagonal elements are negative above 5 GeV, thus violating the 
first inequality (j6]l and providing evidence for breakdown of the independent-particle hypothesis 



^ATLAS 1141 uses a different binning in p, for tfie two particles in tiie pair, thus breaking the symmetry. 

4 



0.25 



15 



0.25 12 3 5 
associated p.j. 









flow fluctuations § 








factorization 






I 








nonflow 1 
































\/ 




















■ 3A 





































































15 0.25 1 2 3 5 15 0.25 1 2 3 5 
associated p.j- associated 



15 



Figure 3: (Color online) Test of inequalities |6j for the correlation matrix V„a(p1 , p'') measured by ALICE 1121 for the 
10% most central Pb-Pb collisions at 2.76 TeV. Green: the second inequality is an equality within errors. Blue: strict 
inequality (flow fluctuations). Red: inequalities violated (congelations cannot be explained by hydrodynamics). From left 
to right: n = 1, 2, 3. Only statistical errors were taken into account. The boundaries of the p, bins are 0.25, 0.5, 0.75, 1, 
1.5, 2, 2.5, 3, 4, 5, 6, 8, 15 (in GeV/c). We use the standard terminology and denote by trigger (associated) the particle 
with the higher (lower) p,, although this distinction is irrelevant here. 



underlying the hydrodynamic picture. The observed effect can be quaUtatively explained by 
correlations from jets: the trigger and associated particles are separated by a pseudorapidity 
gap and cannot belong to the same jet. But they can belong to the same pair of back-to-back 
jets, in which case their azimuthal angles are typically separated by A(p ^ tt, thus giving a 
negative contribution to diagonal elements for odd n. Below 5 GeV, all matrix elements are 
compatible with inequalities (|6]l, except for one point (which is again compatible with a statistical 
fluctuation). As observed for n -2, factorization is broken even at low p,. 

Finally, for n = 1 (Fig. |3] left), the inequalities (|6]l are massively violated, both at low and 
high p,, thus indicating that nonflow effects are important. One of these nonflow effects is the 
correlation induced by global momentum conservation, which only contributes to the first Fourier 
harmonic. After correcting for this effect, data are dominated by collective flow lfT4l[T6l . 

Comparisons between hydro and data so far only address the "single particle v„ip,y\ which 
is inferred from the correlation between a single particle with momentum p, and all the particles 
in a reference detector. Such measurements amount to averaging the correlation matrix. Fig. [3] 
over a line or a column. As we have shown above, the detailed structure of the matrix contains 
much more information; in particular, it provides direct insight into event-by-event fluctuations. 
The correlation matrix defined by Eq. (j5]l can be directly evaluated in event-by-event hydrody- 
namics. Future event-by-event hydrodynamics should address the full structure of two-particle 
correlations, which can be directly compared with experimental data. 

For the sake of brevity, we have only discussed the simplest correlations, i.e., pair corre- 
lations. Much additional information is contained in higher-order correlations. Higher-order 
cumulants give direct information on the magnitude of event-by-event flow fluctuations ifTTl . 
A much wider range of possibilities is opened up by mixed correlations between event planes, 
which have recently been measured at the LHC 11181 [19] and can be directly evaluated in event- 
by-event hydrodynamic calculations ll20ll . Such higher-order correlations are likely to play an 
important part in the near future of hydrodynamic calculations. 
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4. Recent progress and open issues 

It has long been recognized 121 J that the dominant source of uncertainty in comparisons be- 
tween hydro and experimental data lies in the modeling of initial conditions. The modeling of 
event-by-event fluctuations introduces further uncertainty and has triggered much recent activ- 
ity II22I l23l |24| . From the point of view of hydrodynamics, it is important to understand the 
hydrodynamic response to a given initial condition: for instance, to what extent elliptic flow 
represents the response to the initial eccentricity of the interaction region. This hydrodynamic 
response has been studied quantitatively ll25l and the importance of nonlinear terms has been 
pointed out [26]. 

As for the hydrodynamic expansion itself, the dominant source of uncertainty is the value of 
transport coefficients 1271 . Bulk viscosity was recently shown to have a small effect on integrated 
flow observables 1281 l29l . Shear viscosity, on the other hand, is known to reduce anisotropic 
flow: the higher the harmonic, the stronger the effect |30|. But higher harmonics also have 
a larger uncertainty from initial conditions: presently, RHIC data are equally well reproduced 
with a minimal shear viscosity 131] or with ideal hydrodynamics |32|, depending on how initial 
conditions are modeled. Studying ultra-central collisions at the LHC may help reducing the 
uncertainty on the shear viscosity l33l . 

The last stage of the hydrodynamic evolution is the hadronic phase. It is not clear at present 
whether hadronic interactions are strong enough for hydrodynamics to be a valid description of 
this phase. There seems to be a consensus that some hadronic interactions are needed in order to 
match spectra of identified hadrons, but that ideal hydrodynamics fails. A transport calculation 
(hadronic afterburner) reproduces data quite well ll34ll . Alternatively, one can use hydrodynamics 
with a non-vanishing bulk viscosity |35|. 

Table [T] is a list of papers published since the last Quark Matter conference and contain- 
ing numerical hydrodynamic calculations applied to ultrarelativistic nucleus-nucleus collisions. 
The list is by no means exhaustive. In particular, it does not include theoretical developments: 
Progress has been made in understanding the relationship between transport theory and hydro- 
dynamics l50l . which is relevant to our understanding of the hadronic phase. Another important 
development is the relativistic theory of hydrodynamical fluctuations |51 1. The correlations re- 
sulting from such intrinsic fluctuations (421, '52 1, which are typically thought of as "nonflow" 
correlations, could thus eventually be studied within the framework of hydrodynamics l53l . 

5. Perspectives 

Hydrodynamics has been the state-of-the-art approach to the soft sector of nucleus-nucleus 
collisions for at least a decade. It was first used to explain the large value of elliptic flow at RHIC. 
Other azimuthal correlations, such as the "soft ridge" [54 1 , were largely thought to be "nonflow" 
effects, and it took a few years to recognize l2| that they were also naturally described by hydro- 
dynamics. In these proceedings, we have pointed out for the first time that the detailed structure 
of azimuthal correlations 1 12 1 might be explained by hydrodynamics alone below 5 GeV. Hydro- 
dynamics encompasses a wider and wider range of phenomena. 

The field of nucleus-nucleus collisions is characterized by a strong interplay between theory 
and experiment. Most experimental measurements of anisotropic flow are biased by the hydro- 
dynamic picture where flow is a single-particle observable. Elliptic flow is often presented as a 
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Table 1: Recent works (since the Quark Matter 2011 conference, listed in chi'onological order) containing numerical 
hydrodynamic calculations applied to heavy-ion collisions. We indicate by check marks whether or not calculations 
involve event-by-event fluctuations of the initial state; whether they are three dimensional or two-dimensional (with 
Bjorken longitudinal expansion); whether they are ideal or viscous; and finally, whether or not collisions are implemented 
in the hadronic phase (hadronic afterburner). 

"single-particle" observable, even though all measurements are inferred from azimuthal corre- 
lations. Triangular flow was actually discovered [2 | through a thorough analysis of these corre- 
lations. We advocate that experiments measure quantities which are simple and unambiguous. 
The pair correlation ([3]l is a good example of such a simple — yet non-trivial — measurement. 
In turn, theory should directly address what is measured, i.e., correlations. This will allow for a 
more fruitful interplay between theory and experiment. 

Hydrodynamics applied to ultrarelativistic heavy-ion collisions is an active and lively field 
of research, and the recent measurements open up exciting new perspectives. The structure of 
pair correlations will be investigated as a function of transverse momentum and rapidity. The 
recently-measured correlations between event planes ifTSl [T9l wiU also further constrain models 
and tighten our knowledge of the initial state. 
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